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Abstract: We present a phenomenological study of a CP-violating two-Higgs-doublet 
Model with type-II Yukawa couplings at the Large Hadron Collider (LHC). In the light 
of recent LHC data, we focus on the parameter space that survives the current and past 
experimental constraints as well as theoretical bounds on the model. Once the phenomeno- 
logical scenario is set, we analyse the scope of the LHC in exploring this model through 
the discovery of a charged Higgs boson produced in association with a W boson, with the 
former decaying into the lightest neutral Higgs and a second W state, altogether yielding 
a bbW + W~ signature, of which we exploit the Ty + H /_ semileptonic decays. 
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1 Introduction 

One of the main tasks of the Large Hadron Collider (LHC) experiments is to study and 
understand the mechanism of Electroweak Symmetry Breaking (EWSB). Recently, the 
ATLAS and CMS collaborations released the results of the search for the Higgs bosons 
with more than 10 ftr 1 data collected at 7 TeV in 2011 and at 8 TeV in 2012 [1, 2). 
Both experiments have recorded an excess of events above the expected background in 
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different decay channels (mainly 77, ZZ and at a mass near 125 GeV. The excess 

is compatible with the Standard Model (SM) Higgs boson. Complementary evidence is 
also provided by the updated combination of the Higgs searches performed by the CDF 
and DO collaborations at the Tevatron [3]. 

Investigating the Higgs mechanism in the framework of the SM constitutes a major 
effort [4]. However, the minimal choice of Higgs sector is (so far) arbitrary. Even if 
the existence of a scalar resonance, compatible with an SM-like Higgs, has already been 
uncovered by the current data, one must take advantage of the unique opportunity to test 
the phenomenology of more complicated Higgs models. 

Much effort has been dedicated over the years to the study of extended Higgs sectors. In 
this paper we consider the two-Higgs-doublet Model (2HDM) with type-II Yukawa sector. 
This model is one of the most popular extensions of the Higgs sector due to its strong 
connection with a tree-level Minimal Supersymmetric Standard Model (MSSM) [5, 6], 
which is one of the most accredited proposals for solving some theoretical inconsistencies 
of the SM. As is well known, the Higgs sector of the MSSM is quite well constrained in 
terms of the number of free parameters on which the masses depend [7-22] . It is possible 
that the Higgs sector lies in a lower mass range than the superpartners of the SM particles. 
In this regard, the 2HDM should be explored as an effective low-energy MSSM-like Higgs 
sector. 

While a tree-level MSSM Higgs sector is strictly CP-conserving (no mixing is allowed 
between the scalar and pseudo-scalar Higgs components), it has been shown that a CP- 
violating effective Higgs sector could be produced by loop corrections under specific cir- 
cumstances [23] . Accordingly, in this paper we adopt a bottom- up approach by considering 
a CP- violating 2HDM with type-II Yukawa couplings. 

Due to its complicated Yukawa structure, a CP-violating parameter space must be 
carefully constrained by theoretical arguments and experimental data. Therefore, our first 
aim is to provide a detailed analysis of the allowed parameter space in the light of recent 
LHC results. We will show that very little CP-conserving parameter space survives these 
data. This exploration of the allowed parameter space has been addressed recently by 
several authors, from different points of view [24-28]. Where there is overlap, we compare 
our results with those obtained by these authors. 

Regarding phenomenology, the only way to unambiguously probe the existence of a 
Higgs sector with two doublets arises through the discovery of a charged Higgs boson, since 
this particle is the hallmark of such a structure of the Higgs sector. Hence, our second aim 
is to profile a charged Higgs boson in the surviving parameter space via a detailed study 
of its production cross-section and decay Branching Ratios (BRs) . 

Then, our third aim is to study the scope of the LHC in discovering a charged Higgs 
state. In this respect, it is well known that the production of a single charged Higgs state 
at a hadron collider proceeds in association with either top /bottom quarks or scalar / vector 
bosons [29, 30]. By taking into account the recent experimental excess observed by ATLAS 
and CMS [1, 2] and the Tevatron [3] (i.e., corresponding to a light Higgs with a mass of 
« 125 GeV), we propose a search strategy for a charged Higgs boson produced in association 
with a W boson and decaying into a bbW final state. In particular, we show that an 
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appropriate choice of the selection cuts would allow the discovery of such a particle despite 
the considerable ti dominated background. 

This paper is organised as follows: in section 2 we give an overview of the considered 
model, in section 3 we analyse the allowed parameter space in the light of both theoretical 
and experimental constraints, in section 4 we present the main phenomenological results, in 
section 5 we briefly comment on possible future developments, and in section 6 we present 
our conclusions. In the appendix we discuss the decoupling limit for the CP-violating 
type-II 2HDM. 



2 The model 



We describe here our parametrisation of the 2HDM with type-II Yukawa couplings. The 
Higgs sector is defined by the presence of two Higgs doublets, with one ($2) field coupled 
to the u-type quarks, and the other ($1) to the d-type quarks and charged leptons [5]. 
We take the 2HDM potential to be 



V = ^(^i) 2 + ^{^2? + A3(*l*i)(*5*a) 



+ A 4 ($I$ 2 )(^$i) + \ [A 5 ($i$ 2 )" + h.c. 
- j { m ?i(^i $ i)+ [m? 2 ($ t 1 $ 2 )+h.c] +m2 2 (*5* 2 )} • 



it. 



1 r 



(2.1) 



The Z 2 symmetry will be respected by the quartic terms (there are no \q or A7 terms), 
and Flavour- Changing Neutral Currents (FCNCs) are constrained [31]. 
We parametrise the Higgs fields as 
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The real and non-negative Vacuum Expectation Values (VEVs) for the Higgs doublets are 
v\ = vcp and v 2 = vsp, with cp = cos/3 and sp = sin/3, and the ratio defines 



tan (3 



(2.3) 



CP violation is allowed, and it is realised by means of the fact that A5 and ^ 2 2 are 
complex numbers. All three neutral states will then mix, with the physical Higgs particles 
Hi (i = 1, 2, 3) related to the weak fields rjj (j = 1, 2, 3) of Eq. (2.2) by 



(2.4) 



In terms of the non-diagonal mass-squared matrix A4 2 , determined from second derivatives 
of the above potential, we have 
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The 3x3 mixing matrix R governing the neutral sector will be parametrised in terms of 
the angles a%, a 2 and 03 as in [32, 33]: 



R 



( ci c 2 

~(C1 S 2 S3 + SIC3) 
Y -C\ S2 C3 + S1 S3 

where c\ = cosai, s\ = sinai, etc., and 



si c 2 s 2 
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"(Cl S3 + S1 S2 C3) C2C3 
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(2.6) 



(2.7) 



For these angular ranges, we have q > 0, S3 > 0, whereas s\ and S2 may be either positive 
or negative. We will use the terminology "general 2HDM" as a reminder that CP violation 
is allowed. 

With all three masses different, there are three limits of no CP-violation, i.e., with two 
Higgs bosons that are CP-even and one that is odd. In the above notation, the three limits 
are [34]: 



Hx odd 
H 2 odd 
H 3 odd 



a.2 ~ ±7r/2, 0:1,03 arbitrary, 
«2 = 0, 03 = 7r/2, ol\ arbitrary, 
02 = CX3 = 0, a\ arbitrary. 



(2.f 



In the general CP- violating case, the neutral sector is conveniently described by these 
three mixing angles, together with two masses (Mi, M2), tan/3 (the ratio between the two 
Higgs VEVs) and the parameter fi 2 = Rem ? 2 /(2cos/3sin/3). From Eq. (2.5), it follows 
that 

' 2 - (2.9) 



k 



When CP is violated, both (.M 2 )i3 and {Ai 2 ) 2 z will be non-zero. In fact, they are related 
by 

(.M 2 )i3 = tan/3(.M 2 )23. (2.10) 

From these two equations, (2.9) and (2.10), we can determine M3 from Mi, M2, the 
angles (0:1,0:2,0:3) and tan/3 [32]: 

M 2 R 13 (R 12 tan /3- R n ) + M 2 R 23 (R22 tan - R 21 ) 



Mi 



(2.11) 



#33 (#31 -i? 32 tan/3) 

where we impose Mi < M 2 < M3. 

Providing also M H ± and fi 2 , all the A's are consequently determined. Since the left- 
hand side of (2.9) can be expressed in terms of the parameters of the potential (see, for 
example, [35]), we can solve these equations and obtain the A's in terms of the rotation 
matrix, the neutral mass eigenvalues, fi 2 and M H ±. The explicit expressions are given in 
Ref. [34]. 

The interest in allowing for CP violation lies in the fact that it may be helpful for baryo- 
genesis [36]. Also, from a more pragmatic point of view, it opens up a bigger parameter 
space, and allows certain couplings to be larger. 
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2.1 Yukawa and gauge couplings 

For the type-II 2HDM, and for the third generation, the neutral-sector Yukawa couplings 
are (assuming all fields incoming): 

rrjT -igm b 1 . 

2mvK cosp 

Hjtt- } [Rv-i-nctxPRjs]- (2-12) 

z sm p 

Likewise, the charged-Higgs couplings are [5] 

H + bi : — ^ V tb [m b (l + 75) tan/3 + m t {l - 75) cot f3], 

H-fb: — l I V^[m 6 (l-7 5 )tan/3 + m t (l + 7 5 )cot/3]. (2.13) 

2V2 mw 

The HjZZ (HjW + W~) coupling is, relative to that of the SM, given by 

HjZZ {HjW + W-) : [cos pRjx + smPR j2 ], for j = 1. (2.14) 

Note that when H\ is CP-odd (H\ = A), then C2 = and this vector coupling vanishes 
[26]. Finally, the HjH + W~ coupling is given by [35] 

HjH ± W T : |[Ti(sin/3% - cosPR j2 ) + fy 3 ](p£ - pJ)- (2-15) 
3 Constraining the parameter space 

The multi-dimensional type-II 2HDM parameter space is severely restricted by a variety 
of theoretical and experimental constraints, which are discussed in the following. 

3.1 Theoretical constraints 

We impose three classes of theoretical constraints: 

• Positivity: In order to have a stable potential, we impose positivity, V(<&i,<I>2) > 
as I $1 1, 1 $2 1 - ^ 00 [37-39]. Additionally, we must insist on a non-trivial solution to 
Eq. (2.11): M| > and M 2 < M3. While positivity may be satisfied for the given 
parameters of the potential, the considered minimum need not be the global one. 
However, it has been shown that, if a local charge-conserving minimum exists, then 
there can be no charge- breaking minimum [40-42]. Nevertheless, the potential of the 
2HDM can have more than one charge-conserving minimum. We therefore check that 
the minimum obtained is the global one, following the approach of Ref. [43]. 

• Tree-level unitarity: We also impose tree-level unitarity on Higgs-Higgs scattering 
[44-48] . These conditions have a rather dramatic effect at "large" values of tan /3 and 
M H ± , though some tuning of \i can extend the allowed range to larger values of tan /3 
[49]. 
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• Perturbativity: We impose the following upper bound on all A's: 

|Ai| < 4tt£, (3.1) 

with £ = 0.8, meaning |Aj| ^ 10. The effect of this is to restrict large values of the 
masses, unless the soft parameter \i is comparable to Mi and M H ±. 

For illustrations of how these theory constraints cut into the parameter space, see Refs. [34, 
35]. 

3.2 Experimental constraints 

Below, we list the different experimental constraints that are important. The SM predic- 
tions of the flavour observables quoted in this subsection are obtained using Superlso v3.2 
[50, 51]. 

• B — > X s ~f. This rare FCNC inclusive decay receives contributions from the charged 
Higgs boson that can be comparable to the contribution in the SM. Since the 
charged Higgs state always contributes positively to the corresponding BR, it is an 
effective tool to probe the type-II 2HDM. The most up-to-date SM prediction for this 
decay, at the Next-to-Next-to-Leading Order (NNLO), gives [50-65]: 

BR(B X s 7)sm = (3.11 ± 0.22) x 10 -4 , (3.2) 

while the combined experimental value from HFAG points to a larger value [66]: 

BR(B -> X s7 ) exp = (3.55 ± 0.24 ± 0.09) x 10~ 4 . (3.3) 

For type-II Yukawa interactions, which we consider here, light charged Higgs bosons 
are excluded by this observable. The actual limit is sensitive to higher-order QCD 
effects and is of the order of 300 GeV, being more severe at low values of tan j3 
[59, 65, 67, 68]. Recently, a higher-order calculation [69] concludes that the 95% C.L. 
is at 380 GeV. While adopting the more conservative limit of 300 GeV, we shall also 
discuss this more restrictive one. 

• B u — > tv t : In contrast to the b — > 57 transitions, where the charged Higgs state 
participates in loop diagrams, the process B u — > tv t can be mediated by already 
at tree level. Since this decay is helicity suppressed in the SM, whereas there is 
no such suppression for spinless exchange in the limit of large tan/3, these two 
contributions can be of similar magnitude [70] . The 2HDM contribution factorises in 
the ratio Rbtv as compared to the SM value. This decay suffers from uncertainties 
from fs and V u b, and using fg = 192.8 ± 9.9 MeV [71], and the combined value 
I Kb I = (3-92 ± 0.46) x 10~ 3 [72], the SM BR evaluates numerically to [50, 51]: 

BR(B U -> tv t )sm = (1.01 ± 0.29) x 10 -4 . (3.4) 

The SM prediction is compared to the current HFAG value 1 [66] 

BR(B U -»• t^ t ) cxp = (1.64 ± 0.34) x 10~ 4 (3.5) 

lr The latest BaBar result BR(B„ -> rv T ) = (1.83l°;" ±0.24) x 10" 4 [73] is not included in this average. 
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by forming the ratio 



In the framework of the 2HDM this leads to the exclusion of two sectors of the ratio 
tanf3/M H ± [65, 70, 74-78]. 

B — > Dtv t : Compared to B u — > tv t , the semi-leptonic decays B — > Div have the 
advantage of depending on |Vd,|, which is known to greater precision than \V u b\- In 
addition, the BR(i? — > Dtv t ) is about 50 times larger than the BR(i? u — > tu t ) in 
the SM. The experimental determination remains however very complicated due to 
the presence of at least two neutrinos in the final state. The ratio 

_ BR(B -> Dtu t ) 
totr ~ BR ( B ^ D£ue y 

where the 2HDM contributes only to the numerator, allows one to reduce some of 
the theoretical uncertainties. The SM prediction for this ratio is [50, 51] 

^ = (29.7±3)xl0- 2 , (3.8) 

and the most recent experimental result by the BaBar collaboration is [79] 

= (44.0 ± 5.8 ± 4.2) x 10~ 2 . (3.9) 

This ratio is also sensitive to a light charged Higgs boson, and leads to complementary 
constraints to the ones following from B u — > tv t [65, 78, 80-82]. 

D s — > tv t : Constraints on a light charged Higgs can be obtained, competitive with 
those obtained from B u — > tv t [83]. The main uncertainty here is due to the decay 
constant fjj a . The SM prediction for this decay is [50, 51]: 

BR(D S -> Tv T ) SM = (5.11 ± 0.13) x 10~ 2 , (3.10) 

using fjj a = 248 ±2.5 MeV [84], and the current world average of the experimental 
measurements gives [66]: 

BR(D S -> tv t ) cxp = (5.38 ± 0.32) x 10" 2 . (3.11) 



• Bd,s A* A* - : These decays are helicity suppressed in the SM and can receive 
sizeable enhancement or depletion from Higgs-mediated contributions. At large tan/3, 
the non-observation of these decay modes imposes a lower bound on the charged 
Higgs boson mass [85, 86]. The most stringent limits for their BRs were reported 
very recently by the LHCb collaboration [87]: 

BR(B S -> < 4.5 x 10~ 9 , (3.12) 

BR(B d -> < 1.0 x 10~ 9 , (3.13) 
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at 95% C.L. Combining these values with the limits from ATLAS and CMS [88] 
results in even stronger bounds: 



BR(B S ->■ /i + /0 < 4 - 2 x 10 
BR(B d -»■ ii + iT) < 8.1 x 10 



,-9 



(3.14) 
(3.15) 



10 



The SM predictions for these branching ratios are [50, 51, 89]: 



BR(B S -> /T/O = (3.53 ± 0.38) x 10 
BR(B d -»• = (1.1 ± 0.1) x lO -1 ' 



9 



(3.16) 
(3.17) 



with BR(5 S — > fi + [i ) being the more constraining. The largest uncertainty is from 
f Ba , we used f Bs = 234 ± 10 MeV for our evaluation. In the type-II 2HDM, the 
experimental limits can be reached for very large values of the Yukawa couplings and 
small charged Higgs boson masses. The constraining power of B ds —> (J, + (J,~ in this 
study is hence rather limited as compared to the other flavour observables. 

• B° — B° mixing: Due to the possibility of exchange, in addition to W exchange, 
the B° — B° mixing constraint, which is sensitive to the term rrtj cot j3 in the Yukawa 
couplings (2.13), excludes low values of tan j3 and low values of M H ± [65, 90-95]. The 
non-perturbative decay constant fs d and the bag parameter B d which are evaluated 
simultaneously from lattice QCD constitute the largest theoretical uncertainty. 

• Rb'. The branching ratio Rb = T z _^^/Tz-^had would also be affected by Higgs boson 
exchange. The contributions from neutral Higgs bosons to Rb are negligible [35], 
however, charged Higgs boson contributions, via the H^bt Yukawa coupling, as given 
by [96], Eq. (4.2), exclude low values of tan/3 and low M H ±. 

• PP ~^ HjX: Two aspects of the recent neutral Higgs searches at the LHC are 
considered [97, 98]: 

— The production and subsequent decay of a neutral Higgs to 77, around M = 
125 GeV is taken to be within a factor of 2 from the SM rate. Assuming 
the dominant production to be via gluon fusion, this can be approximated as 
0.5 < i? 77 < 2, where we define 



We take into account (1) the modified scalar coupling to the fermion or W in the 
loop, (2) the pseudoscalar contribution, and (3) the charged Higgs contribution 
on the 77 side. This condition (3.18) mainly constrains the Yukawa couplings of 
H±. In particular, the (dominant) Hiti contribution to the loop integrals should 
be comparable to that of the SM, meaning 



F(H X -»• gg)BR(H 1 77) 



(3.18) 



•77 — 



T(H SM -»■ gg)BR(H SM -»■ 77) ' 




in 2 13 tan 2 /3 



P 2 (n) = 0(1), 



(3.19) 
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where P(t) represents the ratio of the pseudoscalar and the scalar contributions 
to the loop integral [6] , with r% = M 2 / (4m. 2 ) . At low tan there is some freedom, 
either sf c| or sf, should be of order unity, whereas at high tan this constraint 
requires ct\ — ±7r/2 and a<i ~ 0. 

The production and subsequent decay, dominantly via ZZ and WW, is con- 
strained in the mass range 130 GeV ^ M ^ 600 GeV. We consider the quantity 

= T(H J ^gg)BR(H J ^ZZ) 
ZZ T(H SM -> gg)BR(H SM ^ ZZ)' 1 ' ] 

for j = 2, 3 and require it to be below the stronger 95% CL obtained by ATLAS 
or CMS. This constraint thus affects the product of the Yukawa and gauge 
couplings of H2 and H% (see Eqs. (2.12) and (2.14)). In the limit of unchanged 
total width, this implies 



f? 2 R 2 



sin 2 tan 2 



[cos0Rji + sm0R j2 ] 2 < 77, j = 2,3, (3.21) 



where rj is the 95% CL on cx/<7sm- For tan/5 of the order of unity, the first factor 
is "small" when is of order unity, whereas the second factor is "small" when 
I Rjs I is of order unity. 

For larger values of tan/3, we may substitute ct\ ~ ±7r/2 and Q2 — from the 
above consideration, whereupon the constraint (3.21) takes the form: 



cos 2 
tan 2 



(s 3 c 3 y < rj, (3.22) 



which is not very strong. In particular, it is automatically satisfied at large 
tan 0. 

For the total widths of H2 and H3 entering in the numerator of (3.20), we include 
also Hj —?• H\Hi and Hj — > HiZ, in addition to the familiar decay modes of 
Hi. The relevant couplings can be found in [99] and [35]. 

• T and S: For the electroweak "precision observables" T and S, we impose the bounds 
|AT| < 0.10, \ AS\ < 0.10 [100], at the 1-a level, within the framework of Refs. [101, 
102]. While S is not very restrictive, T gets a positive contribution from a splitting 
between the masses of charged and neutral Higgs bosons, whereas a pair of neutral 
ones give a negative contribution. 

• Electron Electric Dipole Moment (EDM): The bound on the electron EDM 
constrains the allowed amount of CP violation of the model. We adopt the bound 
[103] (see also [104]): 

|4| ^ 1 x 10" 27 [ecm], (3.23) 

at the 1-a level. The contribution due to neutral Higgs exchange, via the two-loop 
Barr-Zee effect [105], is given by Eq. (3.2) of [104]. 
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In contrast to the MSSM, in the 2HDM, an additional contribution to the muon anoma- 
lous magnetic moment arises only at the two-loop level. Since we are considering heavy 
Higgs bosons (M U M H ± > 100 GeV) therefore, according to [68, 106, 107], the 2HDM 
contribution to the muon anomalous magnetic moment is negligible even for tan (5 as large 
as ~ 40. 

The above constraints are not independent. Therefore, we do not attempt to add their 
contributions to an overall x 2 > but rather require that none of them should be violated by 
more than 2a. The LHC constraints are imposed at the quoted 95% C.L. 

Among these constraints, basically b — > 57 requires M H ± J> 300 (380) GeV, while the 
different S-meson constraints impose additional constraints at low and high values of tan j3 
(the latter are basically excluded anyway, by the unitarity constraints). The T constraint 
prevents the masses of the neutral and charged Higgs bosons from being very different, and 
thus effectively provides a cut-off at high masses. 

3.3 Two scenarios 

In order to develop some intuition for the viable parts of the parameter space, we shall here 
consider two scenarios. In both of them, we take the lightest neutral Higgs boson mass to 
be Mi = 125 GeV. Furthermore: 

• Scenario 1: low tan/3, intermediate different masses M 2 and Mjj± (non- 
decoupling regime): 

tan = 1, 2, 3; Mi = 125 GeV, M 2 = 150, 200, 300, 400, 500 GeV, 

M H ± = 300 - 600 GeV, Li = 200 GeV. (3.24) 

• Scenario 2: high tan /3, heavy degenerate masses A/ 2 and M H ± (decoupling 
regime): 

tan (3 = 5, 10, 20; M x = 125 GeV, M 2 = 400, 600 GeV, 

M H ± ~ M 2 , 11 ~ M 2 . (3.25) 

For scenario 1, we will consider a range of charged Higgs boson masses, from 300 to 
600 GeV (for 700 GeV, only a few viable parameter points are found), and a range of tan j3 
values. We will typically consider small values of tan f3, of the order of 1. High values lead 
(for fixed fi) to conflict with the unitarity constraints. 

Scenario 2 is rather fine-tuned. The masses M 2 and Mjj± have to be very close to li, 
in order to avoid conflict with the unitarity constraints. It is discussed in some detail in 
the appendix. There, it is shown that in addition to the SM-like region of a\ ~ ±7r/2 with 
a 2 ~ (H2 or H3 being odd under CP), there is also another region, with a\ ~ and 
a 2 ~ ±vr/2 (with H x being odd, see [108, 109]). 
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3.4 Studies in a space 
3.4.1 Scenario 1 

We shall first consider separately the LHC constraints on Hi — > 77 (an allowed range) and 
-^2,3 - > W + W~ (an upper bound). For fixed additional input parameters, 

M H ± = 300, 500 GeV, tan/3 = 1, (j, = 200 GeV, (3.26) 

we show in fig. 1 the result of a scan over 10 6 sets of (ai, 02, 0^3) • The blue regions satisfy 
all the constraints discussed above, as well as one of these LHC constraints, whereas the red 
regions do not satisfy the LHC constraint considered, but all other constraints discussed 
in sections 3.1 and 3.2. 




-0.5 ai/7C 0.5 -0.5 ai/7C 0.5 



Figure 1: Allowed regions in the ot\— a.2 parameter space, without (red) and with (blue) 
LHC constraints, for M2 = 300 GeV and the additional parameters given in Eq. (3.26). 
Left, blue: regions surviving the H\ — > 77 constraint. Right, blue: regions surviving the 
H2.3 — > W + W~ constraint. 

The allowed regions are for these parameters rather independent of M H ± in the range 
300-400 GeV, but start shrinking around 500 GeV and vanish around M H ± ~ 600 GeV. (If 
we allow £ = 1, see Eq. (3.1), the allowed values of M H ± reach out to about 700 GeV.) The 
underlying checkered pattern is due to the positivity constraint, together with M| > M| . 

When we impose the LHC constraints discussed above [97, 98] (scanning now over 
5 • 10 6 points), we obtain the blue regions for the effects of the Hi —> 77 (left) and the 
1/2,3 —> W + W~ (right) constraint. For the considered case of tan/3 = 1, M2 = Mjj± = 
300 GeV we see that these only very marginally overlap. 

Imposing then both of these LHC constraints, we show in figure 2 the resulting surviv- 
ing parameter space in green on top of the red regions, now for M2 = 400 GeV, and two 
values of charged-Higgs mass, 300 and 500 GeV. Note that although the LHC experiments 
exclude an SM Higgs with a mass from about 130 GeV to about 600 GeV, there are still 
viable regions of parameter space for the second (and third) Higgs state to be in this region, 
since it may couple more weakly than the SM Higgs boson. 
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Figure 2: Allowed regions in the 01-0:2 parameter space, without (red) and with (green) 
the LHC constraints, for M<i = 400 GeV and the additional parameters given in Eq. (3.26). 



A striking first observation is that the allowed regions are very much reduced, only 
values close to a\ = ±7r / 4 are now allowed for the low value of Mjj± , and a somewhat higher 
value for the higher value of Mjj± . This is a little different from the results reported recently 
for the CP-conserving model, where it was found that only a region around a = is allowed 
[24]. We recall that in the particular CP-conserving limit of 02 = 03 = (corresponding 
to the heaviest one, H3, being odd under CP), a = corresponds to a± = ±ir/2. That 
region is here found to violate unitarity. 




Figure 3: Allowed regions in the 02-03 parameter space, without (red) and with (green) 
the LHC constraints, for the parameters given in Eq. (3.26). 

For a complementary view of the allowed region, we show in fig. 3 the corresponding 
projections onto the 02-0:3 plane. Much of the region with 02 < is populated, whereas 
most of the region with 02 > is excluded. All values of 03 are represented. When we 
impose the LHC constraints, the main characteristic is that 02 becomes more restricted, as 
was also seen in fig. 2, and by Barroso et al [27]. Two regions of a-3-values are represented: 
values close to or it/ 2 (in both cases for 02 close to zero). 

A comparison with Arhrib et al [28] indicates that we find a more constrained region, 
presumably because of our tighter constraint on i? 77 . In addition, we find that some points 
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are excluded because of conflict with the electron EDM. 

For higher values of tan/3, several things change. At some point, also negative values 
of a\ become allowed (not shown), and the allowed ranges of a\ move towards ±7r/2 as 
tan/3 increases. At tan/3 = 2 and 3, the allowed region has shrunk by a factor of 2 to 3, 
compared to tan/3 = 1. 

For lower values of M 2 , the impact of the LHC constraints is more severe, but allowed 
points are found, for example also for M 2 = 150 GeV. For higher values of M 2 , the allowed 
region is restricted to some neighborhood of ai = 7r/4 and a 2 = 0, with 03 close to or 
vr/2. 

In view of the results shown in fig. 2, let us comment on the special limit 

a\ = 7r/4, a 2 small (s 2 — 0). (3.27) 
Then, the rotation matrix can be simplified as 

/ 1 1 \ 

-c 3 c 3 V2s 3 , (3.28) 
V s 3 -S3 V2c 3 J 

and the physical states are related to the "weak" states m as 



R * 

V2 



H 1 = ^=( m + r j2 ), (3.29a) 
H 2 = -j= (-C3771 + c 3 r? 2 + V2s 3 m) , (3.29b) 
H3 = ^= {s3Vi ~ S3V2 + V2c 3 773) . (3.29c) 

We recall that, with type-II Yukawa couplings, rji couples to down-type quarks whereas 
rj2 couples to up-type quarks. Thus, in the limit (3.27), the lightest neutral Higgs boson 
couples coherently to the b- and the i-quarks (with strengths proportional to 1/ cos /3 and 
l/sin/3, respectively, and thus, for tan/3 = 1, with the same strength). The heavier Higgs 
bosons, however, have a CP-even content given by 771 — 772 (note the minus sign). For 
tan/3 = 1 only the pseudoscalar component, 7/3, of these heavier Higgs bosons will couple 
to fermions. We note that two CP-conserving limits are contained in this scenario: H 3 = A 
for a 3 = and H 2 = A for 03 = tt/2. 

3.4.2 Scenario 2 

Again, we start out with an overview of the allowed regions of parameter space in the 
absence of LHC constraints. This is presented in red, in fig. 4, for tan /3 = 10 and two values 
of M<2, namely 400 (left) and 600 (right) GeV. The populated regions are at (01,02) — 
(±7r/2,0) and ~ (0,±7r/2), the decoupling regions which are discussed in the appendix. 
They are seen to shrink considerably as the masses are increased from 400 to 600 GeV. 

When we impose the LHC constraints, the regions near (ai,a 2 ) = (0,±7r/2) are no 
longer allowed. Neither is the region near (a\,a 2 ) ~ (— 7r/2,0). This is a CP-conserving 
limit, also commented on in section 3.4.3. Only a small region near (a±,a 2 ) ~ (7r/2,0) 
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Figure 4: Allowed regions in the ai~«2 parameter space, without (red) and with (green) 
the LHC constraints, for M x = 125 GeV, M 2 = 400 GeV (left), M 2 = 600 GeV (right) and 
M 2 - 10 GeV < M H ± < M 2 + 10 GeV. 



remains. The figure includes a range of values for M H ± within [M 2 — 10 GeV, M 2 + 10 GeV] . 
If we take Mjj± = M 2 , then a 2 has to be slightly different from zero, meaning that CP 
is violated. At the CP-conserving parameter point [a\, a 2 , 03) = (tv/2, 0, 0), unitarity is 
violated. 

3.4.3 CP-conserving limits 

Solutions also exist in CP-conserving limits. With #3 CP-odd (H3 = A), we have per- 
formed scans at tan/3 = 3, 5, 10 and 20, for two heavy-mass cases: (i) H 2 = Mjj± = M = 
400 GeV and (ii) M = 600 GeV. In the absence of the LHC constraint, some range in a\ 
around ±7r/2 (corresponding to a around 0) is populated. Imposing the LHC constraints, 
the range in ot\ is constrained, see fig. 5 for the case tan /3 = 5. This is consistent with 
the allowed regions discussed in section 3.4.1, for a 2 ~ «3 ~ (see figures 2 and 3). For 
increasing values of tan/3 (10 and 20), the allowed bands move towards the edges and 
become more narrow (compare fig. 4). 

When H 2 is CP-odd {H 2 = A) and tan /3 = 5, we find allowed solutions at M = 
400 GeV and 600 GeV. But compared to the case H% = A, the mass range for M3 is more 
constrained. 

We do not find any solution when H\ is CP-odd (Hi = A) . The crucial LHC constraint 
is the quantity i? 77 . Others [25, 26] have argued for an interpretation of the 125 GeV excess 
in terms of a pseudoscalar, but we are not able to confirm this. We find points which satisfy 
all other constraints, but not the LHC constraints. Another recent study agrees with this 
[110]. 
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Figure 5: CP conserving case. Allowed regions in the a\-M^ parameter space, without 
(red) and with (green) the LHC constraints, for Mi = 125 GeV, M 2 = M H ± = 400 GeV 
(left), M 2 = M H ± = 600 GeV (right). 



3.5 Studies in tan [3—M H ± space 
3.5.1 The unitarity constraint 

The unitarity constraint plays an important role in delimiting high values of both tan j3 
and Mjj±. This constraint requires the A's to be small, which to some extent is achieved 
by taking the "soft" mass parameter [i large. In fact, in the co-called decoupling limit, 
discussed for the CP-conserving case in [111], and for the present case in the appendix, one 
can respect the unitarity constraints for large masses, provided \i is tuned to these masses: 

M 2 ~ M 3 ~ M H ± ~ [i. (3.30) 

For moderate values of tan/3 (3 — 5), that limit also requires 

/3 ~ «i, a 2 ~ 0, 03 arbitrary. (3.31) 

For large values of tan /3 5), this evolves into the region (ai, a 2 ) ~ ( 7r /2, 0). Furthermore, 
an additional region opens up for large masses and large tan/3, leading to 

Decoupling 1: (ai, a 2 ) ~ (± 7r /2, 0), (3.32a) 
Decoupling 2: (ai, a 2 ) ~ (0, ±vr/2), (3.32b) 

with 03 arbitrary. Two comments are here in order: (i) because of the periodicity of the 
trigonometric functions, regions at ai ~ — ir/2 and ai ~ +7r/2 are connected; (ii) the SM 
limit requires a\ ~ /?, and is thus contained in the region "Decoupling 1". 

In view of the above discussion, in order to determine the maximally allowed ranges 
of tan (3 and M H ± , we scan over some range in /i, starting at the geometric mean 

fxo = V M 2M H ±. (3.33) 
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Figure 6: Allowed regions in the tan f3-Mjj± parameter space, without (red) and with 
(green) the LHC constraints, for M\ = 125 GeV and four values of M2, as indicated. The 
dashed lines show the recent bound at 380 GeV [69] . 



3.5.2 The experimental constraints 

In fig. 6 we show allowed regions in the tan f3-M H ± plane. Again, the larger red region is 
allowed in the absence of recent LHC results, whereas the green region shows what remains 
compatible with these data. We note some reduction in the range of charged Higgs masses. 
Also, at high tan/3, the masses M2 and M H ± tend to be close, as discussed above. 

The "fractal" appearance of these plots is in part due to the finite number of points 
in the scans. Some could also reflect genuine "islands" in parameter space. 

In fig. 7 we show typical values of M3. Note that for each point in the allowed part of 
this plane, some ranges of a's are allowed (see the previous subsection). Each set of a's 
corresponds to a particular value of M3. The values plotted here are those first encountered 
in a random scan over a's. We see that as M2 and Mjj± increase, also typical values of 
M3 increase. 

By allowing a larger value of the perturbativity cut-off £ of Eq. (3.1), higher masses 
of M H ± would be allowed. For example, £ = 1 permits masses above 600 GeV. Also the 
unitarity and the electroweak parameter T constrain this high-mass region. Which of these 
gives the strongest limits depends on the other parameters. 
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Figure 7: Heaviest neutral Higgs mass, M3 in the tan /3—Mjj± parameter space, with the 
LHC constraints, for M\ = 125 GeV and four values of M%, as indicated. The dashed lines 
show the recent bound at 380 GeV [69]. 



4 Benchmark analysis 

In this section we study the profile of the charged Higgs boson at the LHC in view of 
the allowed parameter space analysis. For this, we start by studying a set of candidate 
benchmark points Pi that allow us to synthesise the main features of the surviving models. 

As we have shown in the previous sections, this model depends on eight parameters. 
However, since the \i parameter does not directly participate in the phenomenology of 
interest, if not specified otherwise, we will consider fi = 200 GeV hereafter. The exception 
will be the high-tan j3 case, where \x has to be carefully tuned in order to respect the 
unitarity constraint. Then, we consider points that pass the constraints, at least in the 
charged Higgs mass range 300 — 600 GeV. We are then left with five parameters: tan /3, 
a±, ot2, «3 and M2. 

In table 1, we list a set of 10 candidate benchmark points: we consider the most 
illustrative four of them for determining cross sections and relevant BRs, the rest will be 
discussed only qualitatively. 

This set of points has distinctive characteristics in the phenomenology, as will be shortly 
made clear. Also, we note that Pg and P\o are nearly CP-conserving, with a<i ~ «3 ~ 0. 
All the others correspond to CP-violating scenarios. 
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Table 1: Benchmark points selected from the allowed parameter space. Masses M 2 and 
allowed range of M H ± are in GeV. For P\-Ps, fi = 200 GeV, whereas for Pg and Pio, 
H = M 2 . 

The model has been implemented through the LanHEP module [112] (see [113] for 
details) and the following analysis has been performed by means of the CalcHEP package 
[114, 115]. Furthermore, we have used the CTEQ6.6M [116] set of five-flavour parton dis- 
tribution functions (PDFs). Due to their relevance at hadron colliders, the effective ggHi, 
77-ffj and jZHi vertices have been implemented by means of a link between CalcHEP and 
LoopTools [117], and the numerical results have been cross-checked against the analytical 
results in [118]. 

4.1 Charged and lightest neutral Higgs bosons: BRs 

This subsection is devoted to an analysis of the BRs of the charged Higgs state in the 
allowed parts of the parameter space. In fact, since we are mainly interested in signatures 
of a charged Higgs boson produced in association with a W boson, which involve model 
dependent couplings, it is of fundamental importance to establish some characteristic fea- 
tures of the BRs for some specific points of parameter space. In this regard, we consider 
four points from table 1 and we determine the most important decay modes. We consider 
only BRs > 10~ 4 , rates below this value are not of phenomenological relevance. Then, 
we have six decay modes: WH\, WH2, WH3, tb, ts, tv t , displayed in fig. 8 for selected 
benchmark points. In addition, we remark that the results are presented for an illustrative 
range of Mjj = 300 — 600 GeV, while recalling that the allowed range is always bounded 
in the range Mjj ~ 380 — 450 GeV for an intermediate choice of tan f3. 

First, we consider two points associated to the choice of tan/3 = 1 in fig. 8, P 2 and 
Pq. The dominant decay mode is always tb, and this feature is even reinforced when the 
masses Mi and M 2 are closer (P2), with respect to the case in which they are well separated 
(Pq). However, it is important to remark that the WHi branching fractions, when allowed 
by phase space, are always ~ O(0. 1) and never smaller than ~ O(0. 01). In particular, if 
M H ± > 400 GeV, then the BR for WH\ is ~ O(0. 1), this assures that the suppression 
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Figure 8: Branching ratios of the charged Higgs vs mass. Four benchmark points have 
been considered, as indicated. 



brought about by this decay mode is never stronger than about an order of magnitude for 
a rather large value of M H ± . 

The result does not hold for the W^H\ case when tan/3 = 2 (see fig. 8, P7 and Pg). 
In fact, it strongly depends on the choice of point in parameter space: for Pg this decay 
mode is suppressed down to ~ O(0. 01), while for P7 its branching franction is restored to 
~ O(O.l) because increasing the mixing between rji and 773 via |a2/*7r| = —0.03 — > —0.07 
increases its CP-odd component, and this effect leads to an enhancement. 

Another feature of the tan/3 = 2 choice is that the W^H2 decay mode is always 
dominant as compared to the tb one when Mjj± is large (though not always allowed), 
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Table 2: BRs of the lightest neutral Higgs, BR(Pi — > X), for a set of benchmark points 
extracted from table 1: Pi, P5, P7, Ps, P97 Pio- 

due to the suppression of the — > tb coupling by a factor ~ 2 plus an always sizeable 
—7- W^H2 coupling. We remark that in this scenario the WH\ BR is anyway ~ 0(0.01) 
or bigger 2 . 

Since we are interested in the phenomenology of the charged Higgs boson produced in 
association to vector bosons, it is important to understand the properties of the lightest 
neutral Higgs. For this, we conclude this subsection by presenting a relevant set of H\ 
BRs. As is clear from table 2, the three most important decay mode are always the 66, gg 
and WW* ones. Since the first is the phenomenologically simplest among the three, we 
will only consider this decay channel for studying combined H\ signatures. 

4.2 Charged Higgs: single production mechanisms 

In this subsection we study single charged Higgs production at the LHC for two choices of 
energy, ^fs = 8 TeV and ^fs = 14 TeV. Considering the partonic amplitudes, we have three 
main mechanisms to produce a single charged Higgs boson from hadrons, i.e., associated 
with bosons (B) or fermions (F): 

B: ( 55 ,66) ^PT^ (fig. 9a); (4.1) 
B: qq'^H ± H i (fig. 9b); (4.2) 

F: gg^H + bi (fig. 9c). (4.3) 

In fig. 9 we show the main partonic contributions to the three production channels. 

Note that the process in fig. 9b will generally be disfavoured for two reasons. First, the 
quark density inside the proton is lower than the gluon density, so this channel is suppressed 
in this respect, owing to the typical masses considered (recall that x 2 oc Mjj ± /s, where 

2 The rv T decay mode, on the other hand, strongly depends on the tan /? value. While it is not a primary 
aim of the present paper to analyse such a scenario, when tan/3 ~ O(10) the BR(H ± — > tv t ) could be 
enhanced up to ~ 0(0.01). 
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Figure 9: Single charged Higgs production channels at parton level. 



^/s = 8 or 14 TeV). Second, the intermediate W boson will be largely off-shell, also inducing 
significant depletion of the production rates. In contrast, the channel in fig. 9a receives 
contributions from both quark and gluon initiated processes and can further be resonant 
in the s-channel (for some of our benchmarks). In principle, the box contribution from 
heavy fermions should be included in the set, for it has been proven in several scenarios 
[119-121] that it can lead to O(10 - 100%) corrections when tan/3 ~ O(10). Still, in 
the phenomenological scenario that we propose to be tested at LHC (tan/3 ~ O(l)) this 
inclusion is totally irrelevant for our conclusions. The mode in fig. 9c also benefits from 
counting on two subchannels, though it is never resonant (as M H ± > mt — rrib for the model 
considered here). 

We consider now the four benchmark points P2, Pq, P7 and Ps from table 1 and in 
fig. 10 we plot the cross sections against the charged Higgs mass for the aforementioned 
production channels. In the first place, we confirm that the associated production with a 
neutral scalar is disfavoured. Secondly, the remaining production mechanisms are always 
within a range of an order of magnitude at most. Again, we remark that the results are 
presented for an illustrative range of Mh = 300 — 600 GeV, while recalling that the allowed 
range is always restricted to the range Mh ~ 380 — 450 GeV for an intermediate choice of 
tan/3. 

As regards the fermion-associated production mechanism of fig. 9c, we remark that it 
only depends on the values of M H ± and tan/3 (see Eq. (2.13)), and there is a considerable 
reduction when moving from tan /3 = 1 to tan (3 = 2 (roughly a factor 2) due to the fact 
that the dominant contribution in the coupling is ~ m^/ tan /3, hence the ratio of VEVs acts 
as a reduction factor. The cross section of the fermion-associated contribution at tan /3 = 1 
is ~ 10 — 10 2 (10 2 — 10 3 ) fb when yfs = 8 (14) TeV and it is mostly inversely proportional 
to tan (3. 
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Figure 10: Cross sections of the single production mechanisms vs the charged Higgs mass 
at y/s = 8 TeV (continuous lines) and y/s = 14 TeV (dashed lines). Four benchmark points 
are considered, as indicated. 



The scope of the fermion-associated production mechanism in extracting a — > Wbb 
signature (see below) has been analysed already in the literature, albeit in the MSSM, see 
[122], and we will revisit it in a CP- violating type-II 2HDM in a future publication. 

Instead, here, we concentrate on vector-boson-associated production. The correspond- 
ing cross sections show a complicated behaviour with respect to different choices of param- 
eters. We start our analysis by considering the channel with a final H^W^ state. From 
fig. 10 (P2 and Pq) we see that a choice of tan /3 = 1 plus a low-Higgs- masses scenario (i-*2 : 
M 2 ~ M 3 ~ 300 - 400 GeV) has a cross section ~ 10 - 10 2 (10 2 - 10 3 ) fb when y/s = 8 (14) 
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TeV, and that it is dominant (competitive) with respect to the fermion-associated produc- 
tion. On the other hand, we see that the choice of a high M2 (Pq) and M3 (> 500 GeV) 
favour the contribution from the parton- level channel gg — > ^=2,3 —> H^W^ 1 proceeding 
through the on-shell -ffi=2,3, and this results in a cross section that is always dominant and 
even enhanced when mw + M H ± < Mj=2,3, i.e. ~ 10 2 — 10 3 (10 3 — 10 4 ) fb when y/s = 8 
(14) TeV. 

These qualitative conclusions also hold when tan (3 = 2. Despite an overall suppression 
by one order of magnitude due to the increased value of tan (3, from fig. 10 (P7) we see 
that on-shell production is not taking place (the H2 and H3 masses are too light). Hence, 
starting from a cross section of ~ 10 2 (10 3 ) fb when \fs = 8 (14) TeV at the lower M H ± scale 
we find a rate of ~ 10 (10 2 ) fb when \fs = 8 (14) TeV for high values of M H ±. However, 
the on-shell production is important for P$, where we can see an interplay between the 
on-shell H2 and H3 (M3 ~ 460 GeV) production at M H ± ~ 300 GeV being realised by a 
"double-shoulder" shaped line. In fact, when M H ± ~ 320 (380) GeV the H 2 (H 3 ) on-shell 
production is switched off. In this framework, the cross section is ~ 10 (10 2 ) fb when 
y/s = 8 (14) TeV, but it is increased by an order of magnitude when M H ± < 380 GeV. 

The high-tan (3 benchmark points give rather low production cross sections. For Pg 
and Pio we find cross sections roughly half and a tenth, respectively, of that for Pg. Since 
the latter, after the cuts discussed below, does not yield any useful signal, we have not 
explored Pg and Pio any further. 

Finally, we briefly comment on the sub-dominant channel with a VF-mediated H Hi = \ % 
final state. This channel is unlikely to have interesting phenomenological implication at the 
LHC, at least at the early running stage: when \/s = 8 TeV and with integrated luminosity 
£ int = 10 fb" 1 the cross section is typically just above the threshold for producing a few 
events. Since it is not competitive with the other production channels, we will not study 
this mechanism. 

In the next subsection we consider the neutral-Higgs-mediated production mechanism 
for analysing the scope of the LHC in discovering such a state in the allowed parameter 
space. 

4.3 pp^-H^W^: significance analysis 

In this subsection we analyse the significance of single charged Higgs boson production in 
association with gauge bosons for the set of benchmark points in table 1 (except Pg and 
Pio). All figures herein refer to an integrated luminosity of 100 fb . 

Among the different charged Higgs decay modes, we have chosen to study — > 
W Hi. The decay chain with the Hi — > bb intermediate decay is numerically favoured, so 
that we adopt it here, hence the complete H^ decay chain is 

P ± W ± H 1 W ± bb. (4.4) 

Therefore, we are interested in the significance (X ~ S/VB) analysis of a 2b + 2W final 
state produced via a single charged Higgs. The most important background at the LHC 
for this final state is top quark pair production. However, we will show that a systematic 
reduction of this background is possible. 
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Notice that one may worry here about the contribution of the fermionic charged Higgs 
decay chain 

H + ^tb^ W + bb, (4.5) 

as it yields an irreducible final state that is identical to the chosen one that could be defined 
as part of either the signal or the background. Under any circumstances, though, we believe 
that our top-mass veto (see below) will render this contribution negligible, so we omit it 
here 3 . 

A 6-tagging efficiency of ~ 70% has been assumed for each 6-(anti)quark in the final 
state, and a full reconstruction efficiency has been assumed with respect to the W bosons. 
Among the possible decay patterns of the two W bosons, the semileptonic one was chosen, 
i.e., one hadronic and one leptonic decay, allowing for the full reconstruction of the events 
(unlike the fully leptonic decay mode) and a neater environment than the fully hadronic 
decay mode. 

Hence, the overall selected process for the signal is the following: 

PP -> W T H ± -> W T W ± H 1 -> W T W ± bb -»• 2j + 2b + H + MET. (4.6) 

For each benchmark point, 2 • 10 4 unweighted events were produced. Regarding the 
top background, 4.5 • 10 6 unweighted events (with generation cuts) have been simulated 
in CalcHEP. For both signal and background the standard set of CTEQ6.6M [116] PDFs 
with scale Q = yfs were employed. For emulating a real LHC-prototype detector, a Gaus- 
sian smearing was included to take into account the electromagnetic energy resolution of 
0.15/\/^ and the hadronic energy resolution of 0.5 /yE. 

We describe now the overall strategy for the background reduction procedure. A 
first set of cuts includes typical detector kinematic acceptances and standard intermediate 
object reconstruction, such as W — > jj and Hi — > bb (cuts 1-3). Further, a t-(anti)quark 
reconstruction is used as "top veto" (cut 4). Led by the consideration that a b quark pair 
stemming from the Higgs boson is boosted (unlike the almost back-to-back pair from ti), 
we define the last cut of the following set (cut 5): 

1) Kinematics: standard detector cuts 

pj > 15 GeV, \rj e \ < 2.5, 

pJ>20GeV, M<3, (4.7) 

|A%| > 0.5, |A%|>0.5; 

with rj the pseudorapidity and AR = \J (Ar/) 2 + (A(f>) 2 . 

2) light Higgs reconstruction: 

\M(bb) - 125 GeV| < 20 GeV ; (4.8) 

3 We will attempt extracting this particular H topology in the context of the CP-violating 2HDM in a 
separate publication, with the aim of improving upon the MSSM results obtained in [123]. 
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3) hadronic W reconstruction (Wh — > jj): 



\M(jj) -80 GeV| < 20 GeV; 

4) top veto: if AR(h, W h ) < AR(b 2 , W h ), then 

M(b x jj) > 200 GeV , M T (b 2 lv) > 200 GeV , 

otherwise 1 -H- 2; 

5) same-hemisphere b quarks: 



(4.9) 



(4.10) 



P&i Pb 2 



|P&il |P&al 



> 0. 



(4.11) 



In table 3 we show the efficiency of the previous set of cuts against the simulated 
background for the P 2 and P4 points of table 1, for two H^ masses. There is a clear 
correlation between the M^± value and the efficiency of the top veto (the most effective 
cut of this set): the higher the mass, the higher the efficiency. 



Cut 


BG events 


BG Eff. (%) 


P 2 events 


P 2 Eff. (%) 


s/Vb 


1 




Kin. 


310291 


100 


54.2, 59.4 


100, 100 


0.1, 0.1 


2 




H\ rec. 


263629 


85.0 


53.8, 59.0 


99.3, 99.3 


0.1, 0.1 


3 




W rec. 


256745 


97.4 


52.2, 57.4 


97.0, 97.3 


0.1, 0.1 


4 




top veto 


1689 


0.7 


15.5, 29.3 


29.7, 51.0 


0.4, 0.7 


5 




same-side 6's 


708 


41.9 


11.8, 22.8 


76.1, 77.9 


0.4, 0.9 


Cut 


BG events 


BG Eff. (%) 


P4 events 


Pi Eff. (%) 


s/Vb 


1 




Kin. 


310291 


100 


356.2, 166.4 


100, 100 


0.6, 0.3 


2 




Hi rec. 


263629 


85.0 


351.5, 165.7 


98.7, 99.6 


0.7, 0.3 


3 




W rec. 


256745 


97.4 


341.9, 160.4 


97.3, 96.8 


0.7, 0.3 


4 




top veto 


1689 


0.7 


41.6, 70.3 


12.2, 43.8 


1.0, 1.7 


5 




same-side 6's 


708 


41.9 


32.8, 54.7 


78.7, 77.7 


1.2, 2.1 



Table 3: Consecutive efficiencies of the cuts imposed on the top quark background and 
(top) on the P 2 point, and (bottom) on the P4 point (with M H ± = 310, 390 GeV) of 
table 1. 

After these rather generic cuts are imposed, more signal-based selections can improve 
the significance. The main consideration of the following analysis is that the charged 
Higgs mass can equivalently be reconstructed by either the invariant mass of the four jets 
(2& + 2j), M(bbjj), or the transverse mass of the b jets, the lepton and the MET, MtQMv). 
Let's focus on the M(bbjj)-MT(bb£v) plane: for the signal, either of the two variables (if 
not both) will always reconstruct the correct charged Higgs boson mass, thus producing a 
cross-like shape in the plane defined by the two masses. In contrast, the background events 
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Figure 11: M(bbjj) vs. MT(bb£v) after cut 5 for (unweighted) point P2, with M H ± = 
375 GeV (red, bottom- left) and M H ± = 525 GeV (green, top-right). In black is the 
(unweighted) top background. The dashed lines show M\\ m = 450 and 600 GeV. 

accumulate at ~ 2mt, as can be seen in fig. 11 in which we adopt an illustrative choice of 
charged Higgs masses. 

The presence of long tails for the signal towards regions where the top background is 
heavily reduced allows us to introduce two specific (and alternative) cuts: 

"squared cut": C squ = max (M(bbjj), M T (bb£v)) > M hm (4.12) 
"single cut": C sng = M T {bbiv) > M hm . (4.13) 

The single cut of eq. (4.13) is applied only on Mxibbiv) because the reduction of the 
top background is higher than if compared to a similar cut on the M(bbjj) for the same 
numerical value of M\ im . 

To determine which is the better of the two proposed strategies and what is the optimal 
value for Mi; m , we studied the effects of C squ and C sng for several values of M\\ m . Results 
are shown in tables 4 and 5 for the points P2 and P4, respectively. 

Clearly, a higher value for Mn m results in an increase of the significance, the top 
background is reduced more than the signal. It is important to note that for low charged 
Higgs masses, C S q U seems to perform better than the single cut. However, this is strickly 
true for M H ± ~ 310 GeV only: if a further selection is imposed, restricting the evaluation 



- 26 - 



Cut 


11 


P 2 = 310 GeV 


P 2 = 390 GeV 


Events 


Events 


c / 

bj V -D 


Events 


D/V-D 


M Um = 450 GeV 


P 

^sng 
Csqu 


66.6 
161.1 


6.6 

in f; 
1U.0 


0.8 
U.o 


12.2 
on i 


1.5 

1.0 


M Um = 500 GeV 


Csng 
Csqu 


A ^ 9 
40. Z 

118.8 


6.0 
9.7 


0.9 
0.9 


n.i 

18.4 


1.6 
1.7 


M um = 550 GeV 


Csng 
Csqu 


30.3 
91.0 


5.1 

8.5 


0.9 
0.9 


9.9 
16.1 


1.8 
1.7 


M lim = 600 GeV 


Csng 
Csqu 


24.9 
63.1 


4.7 
7.7 


1.0 
1.0 


8.9 
14.3 


1.8 
1.8 



Table 4: Comparison between C squ and C sng vs Afn m for surviving events and signifi- 
cance with respect to the background. 



Cut 


tt 


P 4 = 310 GeV 


P 4 = 390 GeV 


Events 


Events 


s/Vb 


Events 


s/Vb 


M Um = 450 GeV 


Csng 
Csqu 


66.6 
161.1 


14.5 
25.8 


1.8 
2.0 


29.0 
47.3 


3.6 
3.7 


M Um = 500 GeV 


Csng 
Csqu 


45.2 
118.8 


12.7 
22.4 


1.9 
2.1 


26.3 
43.0 


3.9 
3.9 


M lim = 550 GeV 


Csng 
Csqu 


30.3 
91.0 


10.8 
19.8 


2.0 
2.1 


23.4 
37.9 


4.2 
4.0 


M Um = 600 GeV 


Csng 
Csqu 


24.9 
63.1 


10.0 
17.7 


2.0 
2.2 


20.3 
33.1 


4.1 
4.2 



Table 5: Comparison between C sn g and C squ vs M^ m for P4: surviving events and signifi- 
cance with respect to the background. 

of the significance to the peak-region only 

peak cut: \M - M H ± \ < 50 GeV , (4.14) 

the significance obtained by imposing C sn g, when calculated for all the other charged Higgs 
boson mass values, is always higher than the one obtained by imposing C squ - Here, M = 
min (M(bbjj),MT(bb£i')^ when eq. (4.12) is employed, while M = M(bbjj) when eq. (4.13) 
is employed. 

For the following analysis, the value Mn m = 600 GeV has been chosen as well as the 
selection C sn g, this choice provides the best significance and a narrower peak while keeping 
a sufficient number of signal events (> 10). Should the surviving signal events be less than 
10, it would then be advisable to choose instead the squared cut C squ for the higher survival 
probability of the signal events (despite the lower significance and the broader peak). 

The invariant mass distributions for the points P2, P3, -P4, -P5, and P7 are plotted in 
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Figs. 12-16, each for two values of the charged Higgs mass. Table 6 collects the results for 
the S/\fB analysis for all points. 





M H ± = 


310 GeV 


M H ± = 


390 GeV 




Events 


S/VS 


Events 


S/VS 


tt 


24.9 


peak 


11.9 




9.9 




P l 


3.8 


0.8 






peak 


2.6 


0.8 






P 2 


4.7 


1.0 


8.8 


1.8 


peak 


3.3 


1.0 


7.3 


2.3 




11.3 


2.3 


22.0 


4.4 


peak 


7.7 


2.3 


17.2 


5.4 


Pi 


10.0 


2.0 


20.3 


4.1 


peak 


7.8 


2.3 


16.0 


5.1 


Ps 


21.1 


4.2 


30.2 


6.1 


peak 


13.9 


4.1 


25.0 


7.9 


Pe 


14.0 


2.8 






peak 


9.4 


2.8 






Pl 


3.1 


0.6 


7.4 


1.5 


peak 


2.8 


0.8 


7.3 


2.3 


Ps 


1.2 


0.2 






peak 


1.2 


0.4 







Table 6: Surviving events and their significance after the single cut of eq. (4.13) and after 
the peak selection of eq. (4.14), for all points of table 1, except Pg and Pio- 

As regards Pi, Pq and Pg, no choice of allowed M H ± produces any appreciable signal 
after the whole set of cuts, hence we will not discuss them any further. 

From fig. 12 we learn that a choice of tan/3 = 1 (though a rather low value of M2 
disfavours the production cross section) is enough to produce a visible signal, even for 
M H ± = 310 GeV with the signal peak lying over the background. However, it is clear that 
the signal suffers from the selection cuts, and this critical situation is eased up only when 
M H ± = 390 GeV. In fact, for higher allowed M H ± masses, the signal contains a conspicuous 
number of events (~ 10 in the peak bin), and it is always above the background. 

If we consider the points P3, P4 and P5 the same circumstances occur: the higher 
M H ±, the more visible and clear the signal with respect to the background. Moreover, 
from P3 (M2 = 350 GeV) to P5 (M2 = 450 GeV) the production cross section increases, as 
shown in figs. 13-15. 

If we consider the point P7 (which has tan (3 = 2) in fig. 16, we do not note any 
change from the previous considerations: even when tan (3 grows we still have a considerable 
production cross section in the allowed parameter space, and the signal would be observable 
with respect to the background, at least for M H ± = 390 GeV. 



- 28 - 



a 

o 
o 



c 

2 

Hi 




S 
o 



c 

III 




100 150 200 250 300 350 400 450 500 550 600 
M(bbjj) (GeV) 



100 150 200 250 300 350 400 450 500 550 600 
M(bbJJ) (GeV) 



Figure 12: Point Pi- Number of events integrated with C in t = 100 fb at y/s = 14 TeV 
vs M(bbjj) for signal (coloured lines) and i-quark background (black histogram, labeled 
"B"). 
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Figure 13: Point P3. Similar to fig. 12. 



5 Possible future scenarios 

We shall here discuss possible future experimental developments, and consider their impli- 
cations for the model, in particular for the proposed benchmarks. The basic question is of 
course: which experimental efforts are required to exclude the 2HDM altogether? We shall 
below adress a couple of LHC-related aspects of this question. 
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Figure 14: Point P\. Similar to fig. 12. 
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Figure 15: Point P5. Similar to fig. 12. 



5.1 Higher and more constrained rates for gg — > Hi — > 77 

Several authors have recently argued that the LHC experiments point to an overall rate for 
pp — > H — > 77 that is somewhat high compared to the SM prediction. In our parameter 
scans in section 3 we generously allowed the ratio i? 77 of Eq. (3.18) to satisfy 0.5 < i? 77 < 
2.0. We shall here briefly comment on how the parameter space is further constrained when 
we only allow the upper range: 

1.5 < i? 77 < 2.0 (5.1) 
For the case shown in Fig. 2, namely tan /3 = 1 and M2 = 400 GeV (and two values of 
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Figure 16: Point P-j. Similar to fig. 12. 



tanp=1 M 2 =400 GeV M^soOGeV tanp=1 M 2 =400 GeV M^soOGeV 




Figure 17: Superimposed on Fig. 2, we show (in blue) the remaining allowed regions 
in the a.\-a.2 parameter space, imposing higher rates for the 77 channel, Eq. (5.1), for 
M2 = 400 GeV and the additional parameters given in Eq. (3.26). 



Mjj±), we show in Fig. 17 how the LHC-allowed region gets constrained. Such a develop- 
ment can have dramatic consequences for the model: benchmark points Pi, Pq, P7 and Ps 
would be excluded. 

In figure 18 we show the remaining allowed regions in the tan (3-Mjj± plane, which 
exhibits a significant reduction, as compared with figure 6. While the scans have lim- 
ited statistics, there is an indication that the remaining allowed parameter space starts 
fragmenting into disconnected regions. 
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Figure 18: Allowed regions in the tan @—Mjj± parameter space, without (red) and with 
(green) the LHC constraints, for the range in i? 77 given by Eq. (5.1), and four values of 
M2, as indicated. The dashed lines show the recent bound at 380 GeV [69]. 



5.2 Tightened upper bound on gg — > #2,3 — > W + W (and ZZ) 

Presumably, the search for an SM-like Higgs will continue in the mass range from around 
130 GeV and up. To a first (rough) approximation, a Higgs in this mass region is produced 
via gluon fusion, and decays via WW (or ZZ) bosons. Assuming these upper bounds are 
tightened, it is interesting to see how the allowed parameter space behaves. In Fig. 19 we 
show how the allowed parameter space shrinks if we assume that the upper bound on a 
Higgs-like particle, represented by the quantity Rzz of Eq. (3.20), is lowered by a factor 
0.5. Again, we see a rather dramatic impact of such a development. 



6 Conclusions 

For the channel pp — > H^W^ — > W + W~bb, we have established a set of 7 benchmarks 
for the CP- violating 2HDM with type-II Yukawa interactions. These points all have M± = 
125 GeV, low tan f3, they all violate CP, and allow for a range of charged-Higgs masses. A 
set of cuts is proposed, that will reduce the ti background to a tolerable level, and allow 
for the detection of a signal in the WW —> jjlv channel. Some of the proposed benchmark 
points lead to enhanced H^W^ production cross sections due to resonant production via 
H2 or H3 in the s-channel. 
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tanp=1 M 2 =400GeV M^soOGeV tanp=1 M 2 =40OGeV M^soOGeV 




Figure 19: Superimposed on Fig. 2, we show (in blue) the remaining allowed regions in the 
a\-a.2 parameter space, assuming that the upper bound on an SM-like Higgs is tightened 
by a factor 0.5 for masses in the range 130-600 GeV, for M<i = 400 GeV and the additional 
parameters given in Eq. (3.26). 

Most of the proposed points are in the interior of some allowed domain in the a space, 
and thus robust with respect to minor modifications of the experimental constraints. Some 
of the benchmark points are vulnerable to a higher value of i? 77 . However, the points P%, 
P3, P4 and P5 are not endangered. 

It should also be noted that the proposed channel only benefits from favourable pro- 
duction cross sections and branching ratio at low values of tan (3. In this region, the charged 
Higgs mass is constrained to the range ~ 380-470 GeV. 
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A Appendix. The decoupling limit 

We shall here explore the so-called decoupling limit, which has been studied for the CP- 
conserving case in [111], where 

M H ± ~ M 3 ~ M 2 » M x . (A.l) 

We shall see that the large masses will all be of order fi. Note, however, that this discussion 
disregards the constraints of positivity, unitarity, etc., that are discussed in section 3. 

A.l tan /3 of 0(1) 

For definiteness, we substitute 

M 3 = M 2 = M (A.2) 

into the expressions for the A's [34] , require them all to be small, and also neglect terms of 
order M 2 compared to M 2 . The conditions related to the different A's are: 

Ai : (cfsjj + s\)M 2 ~ s 2 ^ 2 , (A.3a) 

A 2 : (s{sl + cl)M 2 ~4/i 2 , (A.3b) 

A 3 : ^ciM 2 ~2M| ± -/i 2 , (A.3c) 

C/3S/3 

A 4 : c\M 2 ~ 2M 2 H± - //, (AM) 

ReA 5 : c|M 2 ~/?, (A.3e) 

ImA 5 : (cpci + spS\)c2S2 ^ 0. (A.3f) 

Adding Eqs. (A.3a), (A. 3b) and (A.3e), we get 

M 2 ~ fj 2 , (A A) 

as anticipated. Substituting back into Eqs. (A.3e) and (A. 3d), we find 

c 2 ~ 1, (A.5) 

and 

Mh± ~ fi 2 ~M 2 . (A.6) 
The last equation, Eq. (A.3f), provides two solutions, either 

cos(/3 — ai) ~ 0, or sin(2a2) — 0. (A. 7) 

We note that the angle «3 does not enter in these asymptotic conditions (A. 3), and that 
they are all satisfied for 

a\ ~ (3, a.2 ~ 0, «3 arbitrary. (A. 8) 
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A. 2 Large tan j3 

The case of large values of tan f3 requires special attention. Because of over-all factors 
which were left out in Eq. (A. 3), the first on them, Eq. (A. 3a) must be satisfied to a much 
higher degree than the others (a factor l/c| is involved). This means that the expression 

c\s\ + sj (A.9) 

must be maximised. This requires 

sinai = dbl, sina2 = 0, (A. 10) 

consistent with Eq. (A. 8) (and with H2 or H3 being odd under CP), or 

sinai = 0, sina2 = ±1. (A. 11) 

The latter solution (which corresponds to H\ being odd under CP) is not contained in 
Eq. (A.8). 
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